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Abstract General formalism for the two-component description of the inelastic lepton- 
nucleon scattering in the diffractive region is proposed. Nonperturbative contribution 
to electromagnetic structure functions of a nucleon is described by the modified gener- 
alized vector dominance model containing special cut-off factors restricting the phase 

volume of the initial gq-pairs of virtual photon's fluctuations. Perturbative QCD 
contribution is described by the phenomenological model suggested (in nonunitarized 
form) by Forshaw, Kerley and Shaw. 

Formulae needed for a numerical calculation of photonuclear cross sections inte- 
grated over are presented. It is argued that in the case of the photonuclear cross 
sections at superhigh energies of leptons (E > lO^GeV^), integrated over Q^, the fol- 
lowing two-component scheme is good enough: the nonperturbative contribution is 
approximated by the old parameterization of Bezrukov and Bugaev, and perturba- 
tive one is described by the model of Forshaw, Kerley and Shaw with parameters 
determined from DESY data. Corresponding results of numerical calculations of the 
perturbative part, for the cases of fiA and tA scattering at superhigh energies, are 
given. 

PACS number(s): 13.60.-r, 12.40. Vv 

1 .Introduction 

The photonuclear interaction of leptons is , by definition, the process of the 
inelastic lepton-nucleon or lepton-nucleus scattering, 

l + N ^ I + {hadrons). (1.1) 

This process goes through a virtual photon exchange. In the particular case 
when the four- momentum transferred from the lepton to hadrons is large {\/ —q^ » 
IGeV) the process (1.1) is called the deep inelastic scattering; we are interested 
in the present work in the diffractive region of the kinematic variables (any 
^ 0, transferred energy v{= —go) is large, such that xbj ~ 9^/2(70 is small). 
This region is especially important for applications in cosmic ray physics and 
neutrino astrophysics. The typical problem requiring a knowledge of the pho- 
tonuclear cross section is the study of muon and tau propagation through thick 
layers of matter [1]. High energy muons produced in collisions of cosmic pri- 
maries in atmosphere (the process is Zcr + ZAtm — > tt /x) are detected and 



studied by modern large underground (and underwater) telescopes. As for taus, 
they can be produced in the Earth by extragalactic high energy r— neutrinos 
(such neutrinos must always present in an extragalactic neutrino flux, due to 
Vfj^ < — > v-j- oscillations). 

The process (1.1) is reduced to the absorption of the virtual photon by the 
nucleon, 

-f* + N ^ {hadrons), (1.2) 

and, by optical theorem, is connected with the Compton scattering of a virtual 
photon, 

7* + TV ^ 7* + TV. (1.3) 

The Compton scattering in the diffractive region is described by the vacuum ex- 
change which, in turn, is modelled in QCD by the exchange of two or more gluons 
in a colour singlet state. It becomes possible because, in laboratory system, the 
interaction region has a large longitudinal size, and the photon developes an 
internal structure due to its coupling to quark fields. In the diffractive region, 
i.e. at small Xbj, the diffractive 7* AT— scattering is dominant in the Compton 
amplitude, as compared to the direct contribution due to a photon's bare com- 
ponent (although, the direct process cannot be neglected absolutely, especially 
at high energies of leptons). 

These general statements (about diffractive, hadron-like, nature of 'y*N— 
interaction and even about the two-gluon exchange) are not enough for a quan- 
titative calculation of the Compton amplitude or 7* A/' total cross section. One 
must know also the photon wave function or, by other words, the spectrum of 
hadron-like constituents of the photon and, besides, the amplitudes of the scat- 
tering of these constituents on a nucleon. In pre-QCD era these problems were 
considered, almost exclusively, in terms of vector dominance modcils. These 
models use a hadronic basis for a description of the spectrum of photon's con- 
stituents. It is well known, however, that they can describe the main features of 
the 7* A'— interaction (e.g., the Bjorken scaling) only by using rather unnatural 
assumptions on hadronic amplitudes (see, e.g., [2]). With advent of QCD it be- 
came clear that the description of the photon wave function in terms of purely 

hadronic constituents is not consistent: the (/g-component interacting with a 

target perturbatively must be added. Such a picture arises due to inherent 
QCD properties (asymptotic freedom, colour transparency). 

The aim of the present work is twofold. In the first part of the work we for- 
mulate the two-component model of the electromagnetic structure functions of 
a nucleon separating, in the general form, the nonperturbative part which is de- 
scribed by a generalized vector dominance model. The difference with previous 
works on the subject [3,4] consists just in the use in our approach the generalized 
vector dominance, without limiting oneself by low-mass vector mesons {p,LU, (j)) 
only. At the second part of the work we consider the perturbative component 
and carry out the concrete calculations of photonuclear cross sections for leptons 
of superhigh energies. 

The paper is organized as follows. In sec. 2 the general formulas for the 
two-component structure functions are derived. The hard (nonperturbative) 
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component of the structure functions is considered, in a framework of the colour 
dipole model, in Sec. 3. In Sec. 4. we show how the existing widely used formula 
[5] for the photonuclear cross section must be modified to take into account 
the contribution of the hard component. Results of numerical calculations and 
conclusions are presented in the last section. 



2. Two- component model of structure functions 

A. Two-step picture of the Compton scattering 

The starting point of our consideration is the expression based on the perturba- 
tive QCD and the two-step picture of the 7*p interaction: the first step is the 

7* conversion, and the second step is an interaction of the qq-pa.iT with 

the target proton. The total j*p interaction cross section ( summed over all 

possible final hadronic states) is given by the contribution of the gg-channel in 
the imaginary part of the Compton forward scattering amplitude ( fig.l): 



In this expression A - is the imaginary part of the ( qq )p forward scattering 
qq+p 

amplitude, Q"^ and ^/s are a virtuality of the photon and a j*p center-of-mass 

energy, V7(r,r') is the light cone wave function of grg'-fluctuations of a virtual 
photon with transverse or longitudinal polarization. This wave function depends 
on quark and antiquark helicities (r, r'), quark mass (m/) and quark momentum 
variables z, k± {z = k+/q+ is the fraction of the photon light cone momentum 
carried by the quark, k± is the transverse momentum of the incoming quark) . In 
the following we will need the expressions for the squares of these wave functions 
summed over the quark helicities: 

T *,^' 0. T ,^ o. el{k\k^[z'^ + [l- zf]+ml) 

(t....Q-)<,..,(tx...O') = ^ ^^^^^^^ 

In these formulas M and M' are invariant masses of the incoming and outcoming 
-pairs, 

„ k\^m\ ,2 k'\ + 1711 ^ ^ 
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B.Two gluon exchange approximation 



We suppose, further, that the forward scattering amphtude in Eq.(2.1) can be 
written in the form 



Such a form of the {qq)p amphtude is suggested by the model of the structure 
functions [6] based on the two-gluon exchange approximation of perturbative 
QCD [7-10]. The first term in Eq.(2.4) corresponds to the diagonal transition 
when = M'^; the delta- function factor in this term arises due to zero 
transverse momentum transfer to the quark ( at high energies the value of z is 
frozen during the scattering process, and the interaction can, in a general case, 

change only the transverse momentum of qq pair's particles ). In the case of 
nondiagonal transitions ( the second term in Eq.(2.4)) the scattering amplitude 



depends on the square of the transverse momentum, I j_=kj_ — k±, transferred 
to the quark. 

In spite of the fact that the expressions (2.1), (2.4) are based on the per- 
turbation theory, they can be used for an approximated calculation of the pho- 
toabsorption cross section and electromagnetic structure functions of a nucleon 
at low Q^, where nonperturbative effects are important. It is assumed in such 

an approach that the perturbative fluctuation of the 7* into a qq pair is the 
dominant process and the subsequent colour singlet exchange is realized by two 

gluons only. The problem of the inclusion of higher Fock states (i.e., qq g ) 
in the photon wave function was discussed in the third work of ref.[6]. As for 

the two-gluon exchange appoximation for a description of the elastic {qq)p - 
amplitude, the nonperturbative effects are approximately taken into account by 
introducing an effective mass of the gluons and an effective quark-gluon cou- 
pling constant. The description of soft hadronic processes via the exchange of 
two effective gluons was suggested many years ago in the works [7,8]. 
Below we will use the simple ansatz for the function B{l'j_, s) at small 



The l± dependence of this type is obtained in the models of refs.[6-10] if, 
in addition to the two-gluon exchange, a constituent quark model of the target 
nucleon is used. More exactly, in this case one has [6] 



z', k±, z; s) = {S(k'i_ -ki.) j dH'^B{l'\,s) - B{ll,s)}d{z' - z), 

(2.4) 




B{ll,s) = 




(2.5) 



B{ll,s) 



(2.6) 



2 ' 
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where V{1±) is the GGNN-vertex function, 



y(Zx) = l-e (2.7) 

fiG is the effective mass of the gluon, < > is the mean square radius of a 
nucleon. If hg is small, there is the region, 

^ ' N 

in which B{l\,s) ~ l/l^, as in Eq.(2.5). 

In purely perturbative QCD (colour dipolc picture of the interaction [6,10- 
12], leading agln^ approximation) one obtains for the B{l\, s) function: 

B{ll,s)'^^f{x',ll)as{ll)^^, (2.9) 

, Q2 ^ M2 



s 

;2 



(2.9a) 



Here, f{x,l\) is the unintegrated gluon distribution in a proton. One can 
conclude, comparing Eqs.(2.5) and (2.9), that the ansatz(2.5) is not justified by 
the perturbative QCD calculation. 

C. Separation of perturbative and nonperturbative parts 

Below, in this section, we use the general expressions (2.1-2.4) for a separation 
of soft (nonperturbative) and hard (perturbative) contributions in the structure 

functions. 

A nature of the {qq)p interaction is connected with a transverse size of the 

qq-paii: only for a small r± this nature is, due to the asymptotic freedom, 
perturbative while the pairs with a large enough r± interact nonperturbatively. 
An average transverse separation between particles of the pair is 

'r±= V±^relativer, (2.10) 

where v± ^relative is an relative velocity in the transverse direction, 

— 1/-^ 4- 

Vj_,relative — 2 E E 



E = Jzk^ + kl+ mj, E' = J{1 - z)k^ + kl+ rnj, (2.11) 



and T is the lifetime of the gg-fluctuation, 

1 2k. 
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where fc-y is the virtual photon's three-momentum. 
If follows from Eqs.(2.10-2.12) that 



r_L^ (2-13) 

We use the following criterion of the nonperturbative nature of the interac- 
tion: the {qq)p interaction is nonperturbative if 

r±> r^o = T^, (2.14) 

where k±o is some parameter. This criterion, together with Eqs.(2.3) (where 
the approximation m/ = is used) and Eq.(2.13), leads to the corresponding 
constraint for the variable z. 



(M2 + Q2)2 



according to which the {qq)p interaction is nonperturbative if the gg-pair is 
asymmetric [13], i.e., £■ <C -E" or £" <C E. In turn, the high degree of asymmetry 
corresponds (in accordance with Eqs.(2.3)), at a given value of M, to a small k±_ 
of pair's particles and to their alignment in the direction of the virtual photon 
(the "aligned jet" conjecture had been suggested in the work [14]). 



D.Hard component 

It is convenient, for considering the structure functions in the perturbative do- 
main, to use the impact parameter representation. The resulting formula for 
<^t,l{Q^, s), which follows from Eqs.(2.1-2.4), is well known [6]: 

^TT,L(Q^s) = J dzj d\Mr±,s)\ ^^'^(7^, z, Q') f. (2.15) 

The squares of the photon wave functions in impact parameter space summed 
over helicities are given by the expressions 

I ^^(r^, Q') f= (^)'{[.^ + (l - zf] KfiQ T^)+m)Kl{Q r^)}, 
I i'^ir^.z^ Q') P= ^ E C-^nzHl - zfQ'KliQ r^)}, (2.16) 
Q'^=Q'^z{l-z) + m). (2.16a) 
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The factor a{r±, s) in Eq.(2.15) is a total cross section of the interaction of 
the qq-p&h of the transverse size r± with a target proton, 

a(rx,s) = J dH^{l - e-'^^^^)B{ll,s). (2.17) 

Using Eq.(2.9) for B{1^, s), one obtains (for small values of r±) 



TT^ 1 1 

^{r±,s) —rla,{-^)x',ffG{x'^ff,-^), (2.18) 



where xG{x,Q^) is the gluon structure function of a proton. Approximately 

one has ^ 

^eff ~ - (2.19) 

and[12] 

<7(rx, s) « !^ria«(i-)[xG(a;, -^)] . (2.20) 

Below we will use, instead of Eq.(2.20), the phenomenological model [15] for 
<7{'''±,s), containing except the colour transparency factor r\_ also the cut-off 
factor at large r± (which is necessary for a separation of the perturbative part 
of the structure functions). In such models one has 

a{r±,s) r^rle ^(ris)^". (2.21) 

Note, that s-variable enters the r.h.s. of Eq.(2.21) in combination with 
only, i.e., in the same manner as in the formula (2.20) obtained in perturbative 
QCD. 

E.Soft component 

For considering the nonperturbative part of (Jt,l{Q^, s) we rewrite Eq.(2.1) in 
a form of the double dispersion relation. For this aim we perform the following 
change of variables: 

j dz j d^k± j d^k'j_ = 

= TT [ dzz^(l-z f f dM^ I dM'^ I '^'^'} (2.22) 

J J J J 2MM'Jl-cos'^6 



cos 



2MM'y/l - cos^c, 



■^{M^ + M'-IY), (2.22a) 



72 

: 2 _ '_L 



lY = —^=—, [M - M'Y <l\; <{M + M'y. (2.22&) 
z{l - z) 
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Using this change of variables and the expressions (2.2) for the squares of the 
photon wave functions we have, instead of Eq.(2.1), 



^ / iJtT^ / rf^PT,L(^)-W^^ (M^M'^s), (2.23) 

J M'^ J S qq+p 

Pt{z) = ^[z^ + {1 - z)% (2.23a) 
Pl{z) = ^4z(l - z). (2.236) 

The functions Pt,l{z) are normahzed densities of qq states in the 2:-space. Imag- 
inary parts of the forward {qq)p scattering ampUtudes A^'}' {M"^ , M'^ , s) are 

Qq+p 

obtained using Eqs. (2.4,2.5): 

-7mA^_ {M\M'\s) = -^Nce}{M^diM^ - M'')[ / dH'^B{l'l, s)]- 

S qq+p IZTT^ J 

-ImA^_ {M^,M'^,s) = -^-ImA^ {M^,M'^,s). (2.24a) 

S qq+p^ ' MM' S qq+p^ ' ' ^ V y 

For convenience we included the factor ^^gr in the definition of a longitudinal 

{qq)p amplitude Al- 

It is important to underline that for the derivation of Eqs. (2. 23-2. 24) the 
expression (2.5) for the function B{l'^,s) was used. The factorization of a z- 
dependence and the simple form of this dependence given by the functions of 
Eqs. (2. 23a, b) became possible just due to the assumption that B{l\,s) ~ Ij^. 
The use of this assumption is leqitimate if we consider the nonperturbative 
component of the structure functions. 

Now, for the separation of the nonperturbative contribution in the double 
dispersion relation (2.23) one must take into account the restriction of the qq 
pair's phase space due to the constraint (2.14a). After such a restriction the 
nonperturbative part of <Jt,l{Q'^, s) can be written in the form 



<tt,l{Q\s) = Y, J 



dM'^ 

■X 



M2 + Q2 • 



-J^^^,^(M^M'^fcio)l/mA^:^ (M^M'^.). (2.25) 

-/W + y S qq+p 



8 



The cut-off functions r]T,LiM'^ , M' , k^^) arc obtained in the Appendix A. 

,2 

±0' 



In the region M^, M'^ » A;^q, they are given by the simple expressions: 



r?^(M^M'^ fcio) - vdM^M'^kU ~ 6j^. (2.26) 



The dispersion relation (2.25) can be rewritten by introducing (somewhat 

artificially) the discrete basis for a description of invariant mass spectra of qq 
pairs: 



'2 jji/r'2 



X 



/ ^^!2^/'T,^(M^M'^fcio)-/mA^!^ (M^M'^s), (2.27) 

J M'^+Q^ S qq+p 

PtAM\ M'\ fcio) = Vt,l{M\M'\ klo)6{M^ - M^)S{M'^ - Ml.). 

(2.28) 

Here, Pt,l{M^,M' , fc5.o) spectral functions giving densities of qq - states 
of definite masses. 



F. Generalized vector dominance approach for the soft com- 
ponent 

For a quantitative calculation of the nonperturbative contribution to the struc- 
ture functions we use the approach based on generalized vector dominance 
(GVD) ideas [16,17]. According to hte GVD approach, a gq-pair's quantum 
state which is modified by strong interactions of quark and antiquark can be 
presented by the sum of hadronic (vector meson) states with fixed masses, 

_free _int. — ^ p 

\qq ,M> ^ \qq , M -M^SiM^ - M^) \ V^, M^, > . (2.29) 

Further, it is assumed that the vector mesons are those observed in e+e~ anni- 
hilation process and, correspondingly, the photon-meson coupling constants 
(i.e., the coefficients in the sum in Eq.(2.29)) are expressed through the vector 
meson's leptonic widths. 

Following the experience of vector dominance models [18,19], we will dis- 
tinguish vector mesons of different polarizations. We assume that the meson's 
polarization correlates with a polarization of the initial virtual photon. 

We can write now the basic double dispersion relation of GVDM [16,17], in 
which the cut variables, M and M', are the masses of the incoming and outgoing 
vector meson states in the quasi-elastic forward amplitudes T-^'^ of the vector 
meson-proton scattering: 
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xpgX"(M^M'^s)^-^iJmT-'-(M^M'^.), (2.30) 

and the spectral function p'f^if is a density of the vector meson states. The 
connection between transverse and longitudinal amplitudes is 

hmT^{M\ M'\ s) = ^^hmT^{M\M'\ s), (2.31) 

where parameter ^ takes into account the possible difference in interactions of 
vector mesons with different polarizations. The normalization of these ampli- 
tudes is given by the relations 

hmT^{Ml Ml s) = ay„p(s), (2.32a) 

-ImT^Ml Ml s) = ^^cjvM- (2-326) 

One must note that the introducing of the factor ^ in Eq.(2.31) is at variance 
(if ^ < 1) with the initial two-step formula (2.1) based on the assumption that 

there is the strict factorization (i.e., independence of a {qq)p amplitude on the 
virtual photon's polarization ). 

The crucial assumption of our GVD approach consists in the following: the 
M^- and M'^- dependencies of the spectral functions p'f\^ and p^.L (Eq.(2.28)) 
are the same, i.e., 

p^y{M\ M'^ s) - pt.l{M\ M'^ fcio). (2.33) 

The physical sense of this equality is evident: the hadronic basis is applicable 
for a description of the spectrum of virtual photon fluctuations if (and only if) 

qq pairs produced in the 7* -^qq transition interact with the target proton 
nonperturbativcly. The separation of the nonperturbative component leads to 

the restriction of the phase space of qq pairs and, in turn, to the corresponding 
constraint on the spectral function, i.e., on the density of the vector meson 

states. 

Inserting Eq.(2.33) into Eq.(2.30) one finally obtains 

(2.34) 

This expression differs from a standard GVDM formula only by the presence of 
the cut-off factors r]T,L{Ml M^,, A:^g), the expressions for which are given above. 
These cut-off factors, and the nonperturbative contribution to the structure 
functions given by Eq.(2.34), depend on the parameter r±o(= l/k±o). 

In contrast with the two-component models developed in the works [3,4] our 
treatment of the nonperturbative part contains no restrictions on the vector 
meson mass spectrum: the masses of vector mesons may be arbitrarily large. 
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At the same time, due to a presence of the eut-off factors, it needs no special 
cancellation mechanism (consisting, e.g., in an introducing of unphysically large 
nondiagonal terms ) for ensuring a convergence of the sums in Eq.(2.34) [2]. 

The further consistent developing of the two-component model outlined here 
(first of all, the concrete realization of the modified GVDM based on Eq.(2.34)) 
is a subject of the next work of the present authors. Below, in Sec. 5, we suggest 
the simplified approach for practical calculations of photonuclear cross sections 
at very high lepton energies (> lO^GeV), integrated over Q^, in which the 
nonperturbative component of the structure functions is appoximated by the 
known formulae based on the standard (nondiagonal) GVDM working rather 
well at lepton energies < lO^GeV. 



3. Hard component 

For numerical caculations of the hard component we use the analytical form for 
the function (T{r±,s) suggested in ref.[15]. Namely, 

<j{r^,s) = {air^ + a2rl + a^rl fe'^ (rls)^'' . (3.1) 

The numerical values of the coefficients ai,a2,a3,r^, Xh were found by authors 
of work [15] by fitting the ZEUS data. The main practical aim of our work 
is the calculation of the cross section of the inelastic lepton (muon and tau) 
scattering at extremely high energies (up to lO^Gey). Correspondingly, we need 
structure functions cr(Q^,s) at very low values of Bjorken variable x = Q'^/s 
and, therefore, need extrapolation of HERA data to these values. The Reggc- 
type s-dependence of a{r±,s) in Eq.(3.1) is definitely incorrect at very high 
energies due to the unitarity constraint. So, we are forced to use additional 
assumptions connected with an unitarization procedure. The simplest way of 
the unitarization is the following (see, e.g., [20]). An amplitude for the elastic 

scattering of a gg-pair in impact parameter space which is connected with the 
cross section a{r±, s) by the relation 

a{r±,s) = 2 J (fbj_Ima{r±,b±,s), (3.2) 

is, by assumption, purely imaginary at high energies and can be expressed 
through the opacity function fl: 

a{r^,b^, s) = i{l - e5f^('-^.''^.«)}. (3.3) 

The 6_L-dependence of the opacity function ^l cannot be determined without a 
concrete model. We assume the form [20] 

n{r^,b±,s) = a{r^,s)S{bl), (3.4) 

where the profile function S{b\) is given by the formula [20] 

Sibl) = ^e-'^, (3.5) 
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^ lOGey^. (3.5a) 
The unitarized cross section is given by the integral 

CT"™*(r_L, s) = 2 j (fb^{l - e-|f^('--L'''-L.«)) = 



{'"1 



where C = 0.577. 
Asymptotically 



Ei{-x) |™~ -e-", (3.7) 



so, in the limit of high energies 

cr"""(r±, s) ~ 2TTR^XHlns, (3.8) 
if the s-dependence of a{r±, s) is ~ s'^" . 

4. Cross section of photonuclear interaction 

The differential cross section of the inelastic lepton-nucleon scattering is ex- 
pressed through the electromagnetic structure functions Wi, W2 of nucleons: 

;{{Q^-2mf)Wi{Q^iy)+ 



+{2E{E-u)-^)W2{Q\iy)}, (4.1) 

where E is an lepton energy in laboratory system and m; is the lepton mass. 
The variable v in Eq.(4.1) is related with the 7*p center-of-mass energy ^/s used 
above, 

.= ^x--i±^— (42) 
2M/^i~ 2Mp 2Mp- ^ ' 

The functions Wi , W2 are connected with the structure functions aT,L by the 
simple relations: 

aT(Q^l/) = ^T4-l(g^^^), 

ctt{Q\ y) + c7l{Q\ y) = ^^^^W2iQ\ u), (4.3) 



K = v- « V. (4.3a) 

For problems of cosmic ray physics and high energy particle physics one 
needs the lepton-nucleon inelastic scattering cross section integrated over Q^. 
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Introducing the new variable v = ^ (which is the fraction of the lepton energy 
transferred to newly produced particles) we obtain 



da a f dQ ^' 2M„vE 



v^{E,v,mi) = ^ j ^ /^r x 

Ql 

12 0^2 



Qluin - T^, QLx = ^MpvE. (4.4a) 



As it follows from Eq.(4.3), in the region of small Xbj 



(TL{Q^u) + aT{Q^i^) ~ iyW2{Q^iy)-^ ~ ^, (4.5) 

therefore the contribution of small (near Qmin) dominates in the integral 
in Eq.(4.4). So, as far as we are interested in Q-^— integrated values (it is just 
the case of cosmic ray physics) we can safely use for a theoretical prediction of 
the structure functions <Jt.l{Q'^ t^) the models working well just in the small 
Q^— region. In particular, the structure functions (Jt,l{Q^ jt^) at small and 
not very large v (i.e., in the region where the perturbative component can be 
neglected) are well described by the nondiagonal GVDM [21-23]. The sums over 
vector meson states in this model (analogous those in r.h.s. of Eq.(2.34)) contain 
no cut-off functions, and, instead, the convergence of these sums is provided 
by large cancellations between the diagonal and nondiagonal transitions. This 
model, being essentially nonperturbative, did not describe properly the region 
of very small xbj values and therefore must be used in combination with models 
describing the hard component of structure functions. 

The very simple and convenient for applications parameterization of numer- 
ical calculations of the structure functions within the framework of the non- 
diagonal GVDM was suggested in the paper [5]. The corresponding formulae 
are given in the Appendix B.The simplicity of these parameterizations for the 
structure functions allows for the appoximate analytical integration in the right- 
hand-side of Eq.(4.4). The result is: 

v^{E,v,mi) I ^ ^Aa^p{uy{0.75G,h{x)[^lnil+-^)--^^^ 

av nonpert 2n Qmin "^1 + Q 



mm 



^2 J2 , ^2 + + 7^)] + 0-25[(>< 



)ln{l 



QLn + QLn ™I Qmin ^2 Q 



, , fi.75GsH{x) ' + 0.25-^;n(l + ^)]}, (4.6) 
x=l---H^. (4.6a) 
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This expression is written for a nuclear target. For the case of a nucleon 
target one must put A= 1, Ggh = 1- The values of , m^,, and the definition 
of the function Gsh{x) are given in the Appendix B. 

The formula (4.6) differs from the corresponding result of ref.[5] only by 

2 2 

several new terms which arc proportional to ratios — r and — r. These new terms 
are negligibly small in the case of the muon projectile, but become noticeable 
in the taon's case, especially in the limit of large v. 

The perturbative part of is obtained by substitution in the integral in 

Eq.(4.4) the structure functions given by Eqs.(2.15,2.16) with {qq)p cross section 
o'(?'j_, s) given by Eq.(3.1) (we use the framework of the colour dipole model in 
the variant elaborated in [15]) or, in imitarized form, by Eq.(3.6). 

If the nonunitarized form of (T{r±,s) is used (supposing that the gluon sat- 
uration effects are small through the whole region of lepton energies considered 
in the present work (i.e., up to lO^GeV)), the energy dependence of is 
approximately factorized: 

v^{E,v,mi)\ c^fiv,mi)E^^. (4.7) 

av QCD 

Here, \h is the parameter entering the expressions (2.21,3.1). Correspondingly, 
in this case there is the power law rise with a lepton energy of the QCD part of 
the energy loss coefficient b{E,mi) defined by the integral 

1 

b(E,mi)=NA f v'^{E,v,mi)dv^ E^" . (4.8) 
J dv 



The relative magnitude of the unitarization effects depends, in our phc- 
nomenological approach, on the several model parameters, first of all on param- 
eter R which is supposed to be of the order of target size (see Eq.(3.5a)). The 
smaller parameter R, the sooner logarithmic regime of the energy dependence 
sets in. 

The sum of two contributions given in Eqs.(4.6) and (4.7) is, by terminology 
of cosmic ray physicists, the cross section of the photonuclear interaction: 

da . da . . . da 

v—{E,v,mi) = v—{E,v,mi)\ +v—{E,v,mi}\ . (4.9) 

dv dv nonpert dv QCD 

The values of parameters in the expression (3.1) for the cross section a{r±,s) 
must be determined from the comparison of our model predictions with experi- 
mental data for small xbj ■ The example of such a comparison for the particular 
value XBj = 0.000032 is shown on fig. 2. It is seen from this figure that our 
GVDM prediction for F^ixBj, Q^) is clearly at variance with data even in the 
region of relatively low values of Q^, if XBj is small enough. At the same time, 
one can see that even at very small XBj there is a Q^-region, in which the hard 
component is negligibly small and GVDM works well. 
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The set of parameter values in Eq.(3.1) found from this comparison and used 
for numerical calculations is the following: 



ai = 0.99, a2 = 0.7GeV, = -6.23Gey^ 



r 



l = 0.219GeV~\ Aff = 0.387. (4.10) 

This set differs from that obtaned in ref.[15] only by the smaller value of the 
parameter r*}^. 

5. Results of calculations 

We performed calculations of the photonuclear cross section and energy loss 
coefficient for muons and tau-leptons. 

As is mentioned in the Introduction, muons and taus of super high energies 
(> lO^Ge^) are of natural origin: they appear as a result of interactions of 
superhigh energy cosmic rays in atmosphere (muons) and as a result of or 
Pe oscillations of extragalactic neutrinos on their way to the Earth with 

their subsequent interaction somewhere near the (underground) detector (taus). 
Therefore, we need, in practice, in the case of superhigh energies of leptons, the 
cross section of the inelastic lepton scattering on a nuclear target rather than 
on a nucleon one. 

In the nonperturbative part of the cross section the nuclear shadowing ef- 
fects are taken into account (formula (4.6)). As for the perturbative piece, the 
problem of the nuclear shadowing requires a special investigation which will be 
performed in a separate work. Qualitatively, the shadowing effects are small if 

the exponential term in the integral in the expression for the {qq)Z interaction 
cross section, 

(T -^Sr±, s) = 2 [ dH{l - exp[-a{r^,s)T(b)/2]}, (5.1) 

is close to 1. In this formula T{b) is an optical thickness of the target nucleus 
at an impact parameter b , 

T(b) = J PA{b,z)dz. (5.2) 
The shadowing effects are small if 

a\r^,s)^«a{r^,s)^, (5.3) 

or 
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Here, (Jmax is the maximum value of the {qq)p interaction cross section, ro — 
1.2fin. If ^ = 22 (it is the value which we use in all numerical calculations in 
the present paper), Eq.(5.4) gives the condition 

0-rnaa;(r_L, s) << 40to6. (5.5) 

The example of the calculation of <T(rj_, s) is shown on fig. 3, and it is seen 
that the condition (5.5) is violated only at very high values of s, s > IQ^GeV'^. 
It corresponds to typical lepton energies E = vjv ^ s/{2Mp ■ 0.1) ^ lO^GeV. 
Besides, the condition (5.5) is unnecessarily strong and can be weakened by 
taking into account that the nuclear shadowing effect in perturbative part of 
CF^ strongly depends on because, due to properties of Bessel functions in 
the expressions for photon wave functions (Eqs. 2.16), one has 

Q < 1, (5.6) 

and, correspondingly, 

' ±char ~ Q2' 

(5.7) 

where r±char is some characteristic value of a transverse size of the qq-psdv in the 
nonperturbative interaction. The integral over determining the photonuclear 
cross section is dominated by values near and, therefore, the role 

of (/(j-pairs with large rj_, for which the cross section (T{r±, s) is relatively small, 
is, due to Eq.(5.7), sufficient. Accordingly, shadowing effects cannot be too large. 

Note that the value of r± at which the cross section a{r±,s) has a maximum 
value does not coincide (even approximately) with a value of the parameter 
in the formula (3.1) (it is due to a large value of the parameter as, leading to a 
fast reveal of the aar^-term). So, the typical pair's transverse size corresponding 
to the perturbative interaction is determined, in the parameterization of ref.[15], 
by the combination of two parameters, 03 and r^. 

The results of our calculation of the perturbative part of are shown on 
figs.4,5(for fiA and tA interactions). These results can be used, together with 
the corresponding nonperturbative parts (given by Eq.(4.6)), for calculations of 
a and r propagation through thick layers of matter and for calculations of 
detection probabilities of superhigh energy leptons passing through large un- 
derground telescopes. For completeness, we present on figs. 6, 7 also the nonper- 
turbative parts of All calculations on figs. 4-7 are done for standard rock, 
A = 22. 

For a rough estimate of the photonuclear range of leptons in matter it is 
useful to know the energy loss coefficients, 

1 dE , Na f da , 

(in a case of the perturbative part, if there is no shadowing, ^ A and 
we return to Eq.(4.8)). On figs. 8, 9 we show the results of calculations of such 
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coefficients for standard rock. It is seen that at energies E ^ 10^ — lO^GeV the 
contribution of the perturbative QCD part of the total photonuclear interaction 
in the lepton photonuclear energy losses is dominant (~ 70%). 

The quantitatively similar results for the photonuclear energy loss coefficients 
b{E,mi) were obtained in ref.[l], authors of which used for the description of 
the structure functions ALLM[25] parameterization (based on the Regge ap- 
proach), and in rcf. [26], for smaller energy interval (E < 10^ GeV) and only for 
muons (authors of [26] used for ctt.l the predictions of CKMT[27] model (also 
of Regge type)). Needless to say, that the theoretical approaches of refs. [25,27] 
are completely different from the one presented here. 
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Appendix A: The cut-off functions 

In the case of transversely polarized virtual photons the criterion (2.14) and 
Eqs.(2.3) lead to the change (in the approximation m/ = 0; Q{x) is the step 
function ) 

1 1 

3 



■ J dz[z^ + (1 - zf] -^IJ dz[z^ + {l-z) 





2x 



xe( ^^^5 ^(i-^))e( — ^ , z{i-z)) = 

^2(1 + ^)2 ^M'2(l + ^)2 

= r?T(M^M'^ fcio). (^1) 
In diagonal approximation one obtains 



z^ = \- \ {A2a) 

2 V4 M2(l + ^)2 



If >> one has 



, , r?T(M2,fcio)-3 ■ (A3) 

M2(l + ^)2' M\l + ^Y 



Having in mind that nonperturbative {qq)p interaction amplitudes are es- 
sentially diagonal (it is known that, in the black disc limit, at least, the cross 
sections of nondiagonal diffractive processes should be small), we can approxi- 
mate f]T{M'^,M'^, fc^p) by the formula 

r,T{M\ M'\ fcio) - ^r?T(M2,fcio)7?T(M'2,fcio). (A4) 
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In the case of longitudinally polarized photons the analogue of Eq.(Al) is 

1 1 

I J ^^-4^(1 -'^-"IJ - ^^^^ M^dt^-)^ - .(1 - .)) X 

^ 

x^( ,.,2ft_\^ - ^(1 - ^)) = M'^ fcio), (^5) 

and, in diagonal approximation , 

Similarly with the previous case we approximate r7L(M^, M'^, A;^q) by the for- 
mula 

1 Appendix B: Parameterization of structure func- 
tions in GVDM 

In the GVDM version of rcf. [5], the electromagnetic structure functions of a 
proton are well parameterized by the following expressions: 



m{=OMGeV\ m^ = 1.8GeV^, ^ = 0.25. (S3) 

In these formulae (T~fp{v) is the total photoproduction cross section for a real 
photon of energy v, which was parametrized by the equation [5] 

cr^p{v) = 114.3 + 1.647/n2(0.0213i^) /ife (54) 

{y is in units of GeV). The parameterization (B1,B2) works rather well in 
the small Q^— region (Q^ < 1 GeV^) and in the region of not very small Xbj 
(lO^'' < XBj < 10~^), and does not contradict with the most recent HI and 
ZEUS data. 

The generalization of Eqs.(Bl,B2) for nuclear targets is the following [5]: 

one must introduce the common factor A and multiply the first term in figure 
brackets in right-hand-sides in Eqs.(Bl,B2) by the shadowing function Gsh{x) 
defined by the relation 

Gsh{x) = ^{^-l + e--{l + x)), {B5) 
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X = 0.00282Ai/V^p(i/). (B5a) 

As a result, one obtains 



4 
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Figure 1: The diagram representing schematically the contribution of qq- 
channel in the imaginary part of the Compton forward scattering amplitude. 
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0.1 1 10 

q2 (GeV^) 

Figure 2: The proton structure function F2, for the fixed value of xbj, as a 
function of Q^. 1: hard component (the set of parameters, given by Eq.(4.10), 
is used); 2: GVDM part (formulae (Bl) and (B2) are used); 3: summary curve. 
Experimental points are taken from ref . [24] . 
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r (fm) 

Figure 3: The total cross section of an interaction of the gg-pair of the trans- 
verse size r±_ with a target proton, a{r±,s). 1: Eq.(3.1); 2: the unitarized cross 
section (Eq.(3.6)). s — lO^GeV'^ in both cases. 
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Figure 4: The perturbative part of the differential photonuclear (fiA) cross 
section for A = 22 (integrated over Q^), as a function of v, the fraction of the 
muon energy transferred to hadrons, for several values of an initial energy of 
the muon. The expressions (4.4), (2.15), (2.16), (3.6), and the set of parameters 
(4.10) are used for the calculation. 
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Figure 6: The nonperturbative part of the differential photonuclear {fJ,A) cross 
section calculated using Eq.(4.6), for A = 22. 
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Figure 7: The same as fig. 6, except the change n ^ t. 



28 




10" 



10 



10 

E (GeV) 



10 



10' 



Figure 8: The energy loss coefBcient for a muon, as a function of its energy, for 
standard rock(j4 = 22). The sohd hne is the total coefficient, the dashed line is 
the nonperturbative part calculated using BB formulae [5]. 
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